INTERNATIONAL JOURNAL OF HYDROGEN ENERGY 38 (2013) 1988-1996 



ELSEVIER 


Available online at www.sciencedirect.com 

SciVerse ScienceDirect 


journal homepage: www.elsevier.com/locate/he 


• 

International Journal of 

HYDROGEN 

ENERGY 

»*-- - ' *»*>* wm(li 


HI in 

V A \WmI jal | A NNHta 


© ■ 

[h<£ 



Twin cylinder alpha Stirling engine combined model and 
prototype redesign 

L.S. ScoUo a ’*, P.E. Valdez b , S.R. Santamarina b , M.R. Chini b , J.H. Baron c 

a CONICET. Grupo de Energia Solar - Inst., CEDIAC - Fac., Ingenieria-Universidad Nacional de Cuyo, Centro Universitario, 

Parque General San Martin, GG 405, 5500 Mendoza, Argentina 

h Grupo de Energia Solar - Inst., CEDIAC - Fac., Ingenieria-Universidad Nacional de Cuyo, Centro Universitario, Parque General 
San Martin, CC 405, 5500 Mendoza, Argentina 

c Fac. Ingenieria-Universidad Nacional de Cuyo, Centro Universitario, Parque General San Martin, CC 405, 5500 Mendoza, Argentina 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 26 August 2011 
Received in revised form 
3 January 2012 
Accepted 31 January 2012 


Keywords: 
Dynamic model 
Stirling engine 
Solar energy 
Simulation 
Prototype 
Design 


This article shows the development of a twin cylinder alpha type Stirling engine model, as 
well as the redesign of a prototype, its manufacture and preliminary results. The 
mechanical model is coupled with a thermodynamic model which allows an analysis and 
comparison of the theoretical and experimental behavior of the machine and a better 
understanding of the whole process. On the other side, a prototype redesign is presented. 
Thermal circuit was improved in order to eliminate losses, reduce dead volumes and 
improve heat transfer to the working fluid. Modifications of the seal system and the pistons 
design are included as well. An equation for the sizing of the flywheel is developed and 
finally the results obtained in tests between the new and the previous design are 
compared. 

Copyright © 2012, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights 

reserved. 


1. Introduction 

1.1. Stirling cycle machines. Design and simulation 

methods 

Stirling cycle machines can be used as engines, generators, 
combined heat and power equipment (CHP), cooling machines 
or heat pumps. In these devices, Stirling cycle has the 
advantages of allowing the use of alternative energy sources 
and avoiding environmentally harmful refrigerants. 

Many simulation models for Stirling machines can be 
found in literature. However, most of them study the 
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thermodynamic problem, without linking it to the mechanic 
one [1—3,9] Moreover, Chin-Hsiang Cheng [4] has recently 
presented a model that takes into account this connection, 
which shows in a complete analysis the performance of the 
device. Saravia et al. [5] also present a model including 
mechanical variables. The present work uses the same global 
approach for the development of the model as the one in 
Cheng [4], but developed for a twin cylinder V-90 alpha type 
engine. This paper presents also a redesign and the first 
results of an engine built with the methodology and design 
tools shown in [6,7], Some mechanical aspects of importance 
and the adopted solutions are mentioned. 


Publications, LLC. Published by Elsevier Ltd. All rights reserved. 
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Nomenclature 


a 

(3 

e 

T 
C0 

Qp 

dt 


Variables and measure units 

Angular acceleration, rad/s 2 
Angle between a rod and its piston, rad 
Angle between the crankshaft and the expansion 
cylinder axis, rad 
Time, s 

Angular velocity, rad/s 
Linear piston acceleration, m/s 2 
Time step, s 

Forces applied in the expansion piston and rod, N 
Gravitational acceleration, 9.81 m/s 2 
Forces applied in the compression piston and 
rod, N 

Moment of inertia, kg m 2 
Rod length, m 
Piston length, m 
Mass, kg 
Power, W 
Pressure, Pa 
Heat exchange ,) 

Radius, m 

Gas constant (Helium), J K/kg 
Torque, N*m 
Temperature, K 
Cold source temperature, K 
Hot source temperature, K 
Volume, m 3 

Linear piston velocity, m/s 
Weight, N 


Wo 

x 

Xt, 


I 

L 

L p 

m 

N 

P 

Q 

R 

Rgas 

T 

t 

Tc 

Th 

V 

Up 

W 


Subscripts 


1 

2 

3 

4 

4- 3 

5 

5- 3 
c 

c2 

ca 

ci 

cig 

cil 

e 

el 

f 

G 

h 

k 

n 

P 

q 

r 

t 

tan 


Work, J 

Cylinder axis direction 
Piston position, m 

Perpendicular cylinder axis direction 

Expansion piston 

Compression piston 

Crankshaft-rod joint 

Compression rod-piston joint 

Compression piston rod 

Expansion rod-piston joint 

Expansion piston rod 

Compression 

Compression cylinder 

Crankcase 

Initial load 

Crankshaft 

Cylinder 

Expansion 

Expansion cylinder 

Friction 

Rod center of mass 

Heater 

Cooler 

Rod axis direction 

Piston 

Load 

Regenerator 

Perpendicular to rod axis direction 
Tangential direction 


1.2. Stirling engine , heat pumps and refrigeration 

The engine mode of operation of the machine is more versa¬ 
tile than other engines, since the heat source is external, and 
offers the advantage of adapting to various fuel types and 
alternative energy sources such as wood pellets, hydrogen 
and biogas for example. 

As heat pumps, these machines are capable of transferring 
heat from a low-temperature source to higher temperature 
one through the input of mechanical work. The commercially 
existing air-conditioning equipments use this principle, but 
the disadvantage of their use is that they need refrigerants 
such as HCFCs (R-22), HFCs and others. Here, Stirling cycle 
heat pumps can have a major advantage since they use air, 
nitrogen or helium as working fluids, all of them without 
prejudice to the environment. In that sense, it is of great 
interest the applicability of the Stirling cycle to geothermal 
heat pumps, which absorb heat from a buried coil and deliver 
it inside a room, achieving a coefficient performance (COP) 
with values between 2 and 5.The applications of Stirling cycle 
as cooling machines are varied and become important espe¬ 
cially at low temperatures (below -20 °C) where the tradi¬ 
tional cycles have low performance. 


2. Model 

2.1. Thermodynamic sub-model 

A simplified isothermal type thermodynamic model is used 
here, mainly as a support for the mechanic model in the next 
section. In future versions of the program, the model will 
include the corresponding losses of heat and pressure. 

To calculate the mass present in each machine part, initial 
load pressure is used: 

yyic = PcrV c2 / (Rgas ' tc ) 

YH-k — Pci'Vfc/ (Rgas ' tk ) 

m r = P cl -V r /(Rgas-t r ) ( 1 ) 

^■h — P ci'Via/ (Rgas ' Li) 
m e = Pci'V e i/ (Rgas ' te) 

From mass conservation, the total mass is the sum of all 
them at each moment: 

P (Vei V el V fe V r V h \ 

m = — ( — + — + — + — + —) 2 

Rgas V tc t e t k t r t h J 

Writing it for P gives, 
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m*R 


gas 


Vc Ve Vfe VrVh 

tc t e tfe t r t h 



as a simplification, uniform pressure is considered, 
P = Pel = Pc2. 

The total volume is calculated using: 


V = V c2 + V el + V r + V k + V h (4) 

Using P, the mechanical model will calculate the torque and 
the resulting angular velocity at each step and in consequence 
the volume variationAV. Then, in case of compression: 

Qj? = — Wo c = P-AV (5) 

And during expansion: 

Q h = Wo e = P-AV (6) 


2.2. Dynamic sub-model 

A schematic diagram of an alpha type engine is illustrated in 
Fig. 1. Pistons have a phase angle of tc/2. The assembly of the 
exchangers is also shown. Rods join the crankshaft at a radius 
R cig and an angle 0 from the expansion cylinder axis. 

For this first version of the model, ideal exchangers are 
supposed, so the heater wall temperature and that of the 
expansion cylinder are the same and equal to Th. The same 
applies to the cooler and the compression cylinder with Tc. 
Two more considerations are taken into account: the rod 
center of gravity is on its half length, and all parts are strictly 
rigid. Power is delivered in the main shaft and will be obtained 
from the forces over pistons and rods. 

The equation for the expansion piston position is: 


x p i = L p + Lcos (fa) + Rdgcos 6 
and for the compression one: 

x p2 = Lp + Lcos(fa) + RcigSin 6 
where 


fa 


arcsin 



(7) 

( 8 ) 

(9) 


Regenerator 



and 


02 


arcsin 


cos e 





The volume of the expansion cylinder is: 

V el = [L(l - cos ft) + R C i g (l - cos 0 )]ttI& (11) 

and the volume in the compression cylinder: 

Vc2 = [L(l - cos ft) + Rag (1 - sin 0)]irR^ (12) 

Differentiating equations (7) and (8) with respect to time, the 
linear velocity of both pistons is calculated, 


u p i = -Lw 5 _ 3 sin jS-L - R ci gwsin 0 

(13) 

u p2 = -Lw 4 _ 3 sin /J 2 + Rdgwcos 0 

(14) 

Where: 


dfa wRdgCOS 6 
dt Lcos /? a 

(15) 

d(3 2 wRdgSin 6 

W 4-3 = = T O 

dt Lcos (3 2 

(16) 


In these equations oj is the angular velocity of the flywheel 
(derivative of 0 with respect to time). 

To obtain accelerations, equations (13) and (14) are 
differentiated: 

a pl = -L<t> 5 _ 3 cos fa - La 5 _ 3 sin fa - R cig w 2 cos 0 - R ci gasin 6 (17) 


a p2 = -L<t> 2 _ 3 cos fa - La 4 _ 3 sin fa - R C igW 2 sin 0 + R ci gacos 6 (18) 

where 


R ri 


<^ 5-3 


cig 


Lcos 


(acos 6 - w 2 sin 6 + ww 5 _ 3 tan /? a cos 6) 



a 4 _ 3 = -—( - asin 6 - cj 2 c os 6 + ww 4 _ 3 tan (3 2 sin 6) (20) 

LCOSP2 

and a is dw/dt. 

Torques and forces over the expansion rod and piston are 
analyzed now, following Fig. 2. 

Two convenient reference system are used, the x-y fixed to 
the piston (translation) and the n-t fixed to the rod (rotation- 
translation). 

Torque and force equilibrium equations from Fig. 2 are: 

£t 3 =-T p1+ ^F Gt + LF 5y sinft + LF 5x cos ft 

Awjsing + ft) =0 (21) 

in n 4 direction, 

^F„ = —F c „ + F 3 „ - F 5y cos ft + F 5x sinft +WjCosg + ft) =0 

( 22 ) 

and in t 4 direction, 

^3Ft = F Gt + F 3t + F 5y sin ft + F 5x cos ft - Wjsing + ft) = 0 


Fig. 1 - Schematic alpha type Stirling engine. 


(23) 
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Fig. 3 — Dynamic diagram of compression piston and rod. 


From which F 5x , F 3n and F 3t are obtainable, other expres¬ 
sions are: 

Tpl — h-3^2-3 

Fpi = ITr 2 el (P el - P ca ) 

Fsy — F pi p\ + YYlpiQ.pi 

Fgyi — 2^^ 5 - 3M 5-3 (24) 

F Gt — 2^ m 5-3^5-3 
Wi — YYls-3-g 
\hlpi — YYipi • g 

From all the equations in (24) it is worth saying that F pl is 
the force that the working gas exerts over the piston, as 
a result from the pressure difference between the cylinder and 
the crankcase obtained from the thermodynamic model. Now, 
tangential force over the crankshaft can be calculated: 

f tan — F 3 nSin(^ 1 + 6) + F 3t cos(/? 1 + 6) (25) 

Torque exerted by the expansion system is: 

Ti=F tan -R cig (26) 

Force and torque diagram for the compression system is 
shown in Fig. 3. Torque expression is: 

T 3 —Tp 2 — —Gq t — LG 4y sin fi 2 + LG 4x cos (S 2 

+ ^W lS in(3^ + ft) =0 (27) 

for yi 2 direction, force balance is 

G n = — Gen + G 3n — G 4y cos (3 2 — G 4x sin (S 2 

+ W 2 cos (3 ? + a) = 0 (28) 

and for t 2 direction, 

£G t = Get + G 3t + G 4y sin (3 2 - G 4x cos (3 2 - W 2 sin(^3^ + = 0 

(29) 

from which F 4x , G 3n and G 3t are calculable. Other equations 
are: 


Tp 2 —14-3^4-3 
F p2 = tzy 2 c 2 (P c2 - P ca ) 

G 4 y = Fp 2 p 2 T YYlp 2 Qp 2 

Ggvi — 2 ^ m 4 - 3 ^ 4_ 3 

Gd = -Lm 4 _ 3 a 4 _ 3 
W 2 = m 4 _ 3 -g 
Wp 2 = rrip2 ■ 9 



The tangential force exerted by the compression piston and 
rod over the crankshaft is: 


Gtan — G 3n sin^— — 0 + @ 2 ^ + G 3t cos^— — 6 + /5 2 ^j 

Compression torque is: 


(31) 


T 2 — Gtan - Rcig 


(32) 


Power loss by friction is calculated from a simplified func¬ 
tion obtained in [6]: 


N f = 3 x 10' • ( 


3,0154 


(33) 


So the friction torque can be calculated as: 


Nf /30\ 3,0154 

Tr = A = 3 X 10 7 - (—] -w 2t 

CO \7T J 


(34) 


The crankshaft torque is: 


Tcig=Icig-a (35) 

If I is the moment of inertia of the flywheel, then 

^ = Tl + T 2+ T dg -T ; -T q (36) 

Where T q is external load torque, given by the connected 
mechanism. From this expression oj can be calculated at any 
time step. 
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Model flow chart of the computational process is shown in 
Fig. 4. Where: 

w i+1 = J + i- (Ta + T 2 + T c j g -Tf- Tg) (37) 

6 l+1 = e [ + w i+1 di (38) 

a I+1 = (w I+1 — w 1 ) / dr (39) 


Table 1 - Specification of engine main parameters. 

Main specifications of the engine 

Power 

428.22 

W 

Reference pressure 

0.66 

Mpa 

Swept volume 

439.707 

cm 3 

Frequency 

13.83 (830) 

1/s (rpm) 

Beale number 

0.11 


Phase angle 

90° 


Working gas 

Helium 



3. Prototype redesign and fabrication process 

A previously built prototype published in [8] was redesigned 
and constructed in order to improve its performance. Several 
changes were made including the external heat supply device 
for tests, heater design, cooler design, regenerator, flywheel 
size, piston design, static and dynamic sealing mechanisms. 
Main specifications of the engine are shown in Table 1. 

3.1. External heat supply for tests 

First an electric oven was designed and calculated according 
to the power required to achieve the design temperatures of 


THERMODYNAMIC 

MODEL 


/-\ 

START 

V_> 


Initial 

values 


r ‘ +1 = r’ + dr 





V V 

v c2> y el 



/* 




END 






the engine by radiation (Fig. 5). A PID controller was installed 
for the manipulation of the power supplied, controlling the 
gas temperature inside the expansion exchanger. To achieve 
this, PT100 temperature sensors were incorporated inside the 
tubes of the expansion exchanger. 

After performing the preliminary tests, it was noticed that 
the heat transfer was not the expected one and that there 
were dead spots of radiation. Therefore a gas burner was 
designed with the aim of achieving heat transfer mainly by 
convection. This alternative managed to maintain a more 
uniform temperature throughout the expansion exchanger. 
See Fig. 6. 

3.2. Heater , cooler and regenerator sizing 

Details of the design process which was made by a scaling 
methodology can be seen in [6]. The resulting specification of 
the exchangers is shown in Table 2. 

3.3. Improvements to the heater 

After construction and testing it was noted that there was 
a deficit of heat exchange area in the outside of expansion 
exchanger. To solve this problem fins were added to increase 
the exchange area as shown in Fig. 7. Fig. 6 shows the burner- 
expansion exchanger set with fins in operation. 

3.4. Flywheel sizing 

An important instance in improvements was to adopt a 
flywheel of greater inertia than the one used in [8], since it was 



Fig. 4 - Model flow chart. 


Fig. 5 - Electric oven used as external heat source in tests. 
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Fig. 6 — Expansion exchanger with burner incorporated 
inside the combustion chamber. 



Fig. 7 — Engine with new finned expansion exchanger. 


observed that the existing one did not store enough energy to 
overcome the work of compression and mechanical losses of 
a complete machine cycle. 

The calculation described below verifies this. From this 
result, it was proceeded to find a flywheel on the market, 
whose moment of inertia was equal to or greater than the 
estimated one. After this, a new shaft design for the coupling 
of the engine to the electric starter was necessary. This 
included the selection of bearings, shaft coupling and the 
flywheel support design. 

The flywheel of kinematic Stirling engines plays a very 
important role, which is to perform the compression work 
required and also to overcome the friction work loss in that 
part of the cycle where the machine is not delivering power. In 
many engines this is not taken into account and on certain 
occasions it can cause the engine to fail to keep functioning 
despite having a good thermodynamic performance. That is 
why in this work a specific equation to calculate the required 
moment of inertia of the flywheel is developed. The 
compression work, assuming it isothermal is: 


Table 2 — Specification of engine exchangers. 

Compression exchanger 


Number of tubes 

43 

Length of tubes 

181 mm 

Int. diameter of tubes 

3.45 mm 

Expansion exchanger 


Number of tubes 

10 

Length of tubes 

446.2 mm 

Int. diameter of tubes 

4.35 mm 

Regenerator 


Length 

54.26 mm 

Hydraulic radius 

0.117 mm 

Free flow area 

2474 mm 2 

Wire diameter of mesh 

0.21 mm 

Mesh number 

1.874 

Number of layers 

128 

Mesh surface 

3558 mm 2 

Mesh layers diameter 

67.57 mm 


V2 



VI 

Considering the gas as ideal, P 


W c = 


V2 

j MV = mRT 



VI 


VI 


mRT 

V 




W c = mRTln(V 2 /V 1 ) (42) 

In this case, taking into account a charge pressure of 6 bar 
the parameters are calculated in Table 3. 

The kinetic energy K stored in the flywheel is: 

K = ^lio 2 (43) 

And, 


W c = -AK 

So the moment of inertia can be calculated as: 



2AK _ 2W C 

(Ci>2 — (*)\) (Ci>2 — Ct>i) 



The required moment of inertia will be the one that allows 
an acceptable range of angular velocity variation. If the 
acceptable relative angular velocity variation is defined as: 


Aw r 


Aw Ct>2 — 

(x) (O i 



Table 3 
work. 

- Gas characteristics and resulting compression 

m 

7,82E-04 kg 

R 

2078.6 J/kg K 

T 

300 K 

v 2 

8.1325e-004 m 3 

Vi 

3.8021e-004 m 3 

w c 

-371,00J 
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Wherew is the desired average angular velocity in rad/s. Now 
o) 2 can be replaced in equation (45)with its expression from 
equation (46). Operating and making co = o> 1 the equation for 
the required moment of inertia of the flywheel is: 


—2W C 

o) 2 (A o) 2 + 2Aw r ) 



With the engine data and considering the value of the 
necessary moment of inertia is 0.91 kg.m2. 

So far the friction work W f was not considered. It can be 
calculated from the friction power Nf as: 


The calculation of the moment of inertia including the 
friction work is then: 


-2(w c+ ^ 
r _ V h 

oj 2 (Aw^ + 2Aw r ) 



In this case the calculated friction power is 200 W, so the 
needed moment of inertia is 0.93 kg.m2. It is thus noted that the 
friction work does not significantly influence the calculation. 


3.5. Redesign of pistons and seals 


N f 

W f = —dd (48) 

0) 

Where the value of to consider is radians, considering 
compression occupies about half of the cycle approximately. 
Then, 



rcNf 

0 ) 



The pistons of the previous engine [8] were the original ones 
of the base machine used (an air compressor). Improvement 
in the design of the expansion piston was decided in order 
to increase the distance between the heat source and 
piston rings, since high temperatures can affect the sealing 
materials. 

In Fig. 8 the stainless steel extension drawing of expansion 
piston can be seen. Inside it, horizontal baffles were placed in 


in 
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>o 

o' 

m 
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in 

o 

ii 

a j 
>- 


Fig. 8 - Expansion piston. Design detail. 


































































INTERNATIONAL JOURNAL OL HYDROGEN ENERGY 38 (2013) 1988-1996 


1995 



Fig. 9 - Reinforced PTFE rings. Fig- H - Evolution graphic comparing the temperature of 

_ the expansion exchanger tubes in the new and old design 

respectively. 


order to prevent the heat convection and radiation between 
the hot part (top) and the cold one (bottom). 

As regards the piston seal, in the earlier constructed device 
[8] it had been assembled with the original rings of the air 
compressor that are made of steel. They allowed a great 
amount of gas passing from the cylinders to the crankcase 
during compression and expansion. 

To reduce such losses of gas to the crankcase, reinforced 
PTFE rings backed up by inside o-rings were designed. Fig. 9 
shows the ring that works as a dynamic seal. Likewise, 
a PTFE guide ring was added for the piston to move smoother 
and avoid erosion. Fig. 10 shows the overall assembly of the 
expansion piston of the thermal machine. 


3.6. Tests 

Tests were performed for the equipment working as a refrig¬ 
erator in order to make it better comparable with previous 
tests. During the test the temperature evolution in the 
expansion exchanger’s wall was measured as the engine was 
moved by an auxiliary electric motor. 

This test was conducted for 9 min. The following graph 
shows this evolution compared with that obtained in [8] for 



Fig. 10 - View in detail of the piston rings assembled. 


the previous design of the engine. It should be mentioned that 
they are not directly comparable because the one from [2] was 
performed at 7 bar and 200rpm while the current one was held 
at 2.4 bar and 830 rpm. Furthermore, 1L of water was placed in 
a flexible container in contact with the exchanger’s tubes. It 
can be seen in the graph of Fig. 11 that the temperature 
reached -23 °C. 


4. Results 

As regards the model developed, Fig. 12 shows angular 
velocity variation in time, for different flywheel moments of 
inertia (6 kg-m2; 0.6 kg-m2; 0.06 kg-m2), all of them for an 
initial angular velocity of 850 rpm. If less inertia is used, the 
angular velocity decreases rapidly, meaning that the engine 
does not work for this pair of I and cj. 

For greater moments of inertia, angular velocity estab¬ 
lishes in a longer time, but its amplitude is lower, representing 
a better running. As the steady-state regime is reached, the 
average angular velocity approaches the same value, which is 
near 1900rpm, for all values of I used in Fig. 12. 


Evolution of angular velocity with time 
2500 


2000 


E 

e- 


& 1500 

o 

_o 

a> 

> 


TO 

a> 

c 

< 


1000 



6 8 10 

Time (s) 



Fig. 12 — Angular velocity as a function of time, for different 
flywheel moments of inertia. 
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Fig. 13 - Angular velocity in time, for different initial 
velocities. 


Fig. 13 shows angular velocity variation in time, for 
a flywheel moment of inertia of 0.06 kg* m2, but different 
initial velocities (1500rpm and 2500rpm). 

For initial rotational velocities under 850rpm, velocity 
function decreases fast. Above that initial condition, regard¬ 
less of its value, oj establishes on 1900 rpm. Amplitudes, after 
transient, are both equal. 


5. Conclusions 

The elaboration process of the model allows a better under¬ 
standing of the system as a whole. The thermodynamic model 
used is a basic isothermal one, with results compatible with 
those of Schmidt analysis. Its simplicity is used to determine if 
the mechanical model is running as it should. Once this is 
checked, an improvement of thermodynamic model should be 
done, in order to approach its results to the experimental 
ones, and make more reliable predictions. 

Next version of this model should solve the influence of 
variables such as piston, rod and crankshaft weight over 
performance, rotational velocity and power, in connection 
with pressure loss, heat loss and regenerator performance. 

With reference to the prototype, as mentioned above, the 
choice of a gas burner as the heat source initially gives good 
results, but it is to consider the possibility of increasing the 
airflow to improve heat transfer by convection to the finned 


tubes. The addition of fins to the expansion exchanger tubes 
increased the heat transfer area, producing good results, as it 
turned out that the temperature was more distributed in the 
tubes as shown in Fig. 6 and also an increase in the rate of heat 
transfer. 

The developed equation for the flywheel allows a quick and 
simple sizing for future applications. 

The redesign work on pistons and seals showed very good 
results, as it was possible to measure a pressure difference 
between the thermodynamic circuit and the crankcase, which 
demonstrates that the new reinforced PTFE rings are working 
as expected. Measurements of mechanical friction should be 
carried out. As for static seals, their performances were very 
satisfactory since working fluid losses were not detected. 
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